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Abstract 



We analyze the peculiarities of the motional narrowing effect in disor- 
dered one-dimensional Frenkel chains with off-diagonal disorder, induced by 
uncorrelated Gaussian fluctuations in the positions of the host units. A clear 
difference in the scaling laws with respect to the magnitude of the positional 
disorder and the chain size is found, as compared to those for uncorrelated 
diagonal disorder as well as for off-diagonal disorder modeled by uncorrelated 
randomness in the nearest-neighbor couplings. The origin of such a difference 
is discussed in detail. 

PACS number(s): 71.35.Aa; 36.20.Kd; 78.30.Ly 
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I. INTRODUCTION 



The concept of motional narrowing, earlier raised by Knapp [|T| for the on-site energy (or 
diagonal) disorder, has been very fruitful in explaining many optical phenomena in quasi- 
one-dimensional systems like J aggregates of polymethine dyes and conjugated polymers 
(for a review, see Refs. 0-^] and references therein). This effect is manifested as a decrease 
of the magnitude of the diagonal disorder as soon as the states of individual molecules are 
collectivized due to the intermolecular interaction and form excitonic states. The decrease 
depends on whether or not the disorder is small enough to be regarded as a perturbation. In 
the perturbative case, the suppression factor is determined by the square root of the whole 
number of molecules in an aggregate, while in the nonperturbative case this number should 
be substituted by the number of coherently bound molecules A similar narrowing effect, 
but different suppression factor, was recently reported for the case of dynamic disorder ^-0. 

It was found in Ref. that the numerically simulated absorption spectrum of polysilane 
with an uncorrelated Gaussian distribution of nearest-neighbor couplings are similar to those 
for an uncorrelated diagonal disorder (see also 0). In contrast, numerical simulations of 
off-diagonal disorder given by Gaussian randomness in the molecular positions PJl^l found 
that the different behavior of the optical observables was very different from that expected 
from standard motional narrowing arguments |T[]. 

The main goal of the present Letter is to uncover the origin of such a difference. It will 
be shown that, in the case of configurational disorder, certain correlations appear in the 
distribution of hopping integrals, in spite of the fact that the distribution of the molecular 
positions is uncorrelated. This finally results in a scaling law for the typical fluctuation 
of the Frenkel Hamiltonian matrix with respect to the magnitude of positional disorder 
and the number of molecules different from that for uncorrelated diagonal disorder. As a 
consequence, the main features of the exciton optical response are largely affected. 
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II. DESCRIPTION OF THE MODEL 



Consider a collection of two-level molecules forming a one-dimensional lattice. For our 
present purposes, we will neglect the static inhomogeneous offset energy of the molecules 
imposed by the surrounding host medium (diagonal disorder). Under this assumption, the 
effective Frenkel Hamiltonian describing the system can be written in the form 

^ = II Jnm\n){m\ . (1) 

nm 

Here, the state vector \n) denotes the n-th molecule being excited and the site index n 
lies within the symmetric domain — (A^ — l)/2 < n < (A^ — l)/2 (A^ assumed to be odd). 
Jnm is the intersite interaction, which will be considered to be of dipole origin, so that 
Jnm = J /\n ~ m + in ~ 'Cm|^5 fi ^ m, and Jnn = 0; ^„ is the deviation of the n-th. molecule 
from its regular position, which we suppose that it occurs only along the directions toward 
the two adjacent molecules. This restriction allows us to replace the fluctuations of vector 
positions by scalars. We note that, arguably, this assumption does not affect the conclusions 
and is introduced only to simplify the analytical treatment. The distribution function of 
is chosen to be Gaussian, 

nin) = A^e^v{-$^\ , (2) 



with variance o"^. 

To study the motional narrowing effect, one should rewrite the Hamiltonian (|l]) in the 
excitonic representation using the eigenfunctions of the unperturbed Hamiltonian (with no 
disorder). For the sake of simplicity, we assume periodic boundary conditions. Then, Bloch 
plane waves are the proper eigenfunctions of (|l]) in the absence of disorder: 

\K) = ^ (3) 



'N n 

where K = 2nk/N belongs to the first Brillouin zone (-(A^ - l)/2 < A; < (A^ - l)/2). 
In the i^T-representation, the Hamiltonian (|1]) reads 
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'H = Y.Ek\K){K\+ ^kk'\K){K'\ , (4a) 

K K,K' 

EK = 2jy 4 cosiKn) , (4b) 

^KK' = ^ E 5 Jn^e^^^"-'"''"^ , (4c) 
SJnm = j(-. ^ ' (4d) 

where the scattering matrix A^k' is real due to the symmetry of the summation region in 
(^. Its diagonal part gives rise to exciton energy shift, yielding inhomogeneous broadening 
of exciton levels for an ensemble of chains. This effect appears to be the main consequence 
of the presence of disorder if the off-diagonal term can be regarded as a perturbation. The 
typical fluctuation of A^k has a direct relation with the inhomogeneous width of the cor- 
responding exciton state \K). 

For nonperturbative magnitudes of the disorder, the off-diagonal elements A^k' {K 7^ 
K') mix the exciton states, resulting in their localization on chain segments of a typical 
size smaller than the chain length and subsequently affecting the exciton optical response. 
Recall that, for a perfect chain, only the state with i^' = is coupled to the light and carries 
the entire exciton oscillator strength, which is then times larger than that for an isolated 
molecule. Being mixed with other (nonradiative) states {K 7^ 0), the radiative state loses 
a part of the oscillator strength due to its spreading over to the nonradiative ones. Thus, 
an effective number, usually called number of coherently bound molecules A^* < A^, replaces 
the system size as the enhancement factor of the oscillator strength of the localized exciton 
states [0]. It reflects the typical number of sites on which the localized exciton wave functions 
have a signiflcant magnitude or, in other words, the number of molecules within a typical 
localization segment. Accordingly, the inhomogeneous width of the optical exciton line will 
also be subjected to renormalization [0]. In the next Section, we will formulate a condition 
separating the perturbative and nonperturbative ranges of disorder magnitudes in order to 
treat the exciton optical response. 
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III. MOTIONAL NARROWING EFFECT 



To gain insight into the magnitude of the typical fluctuation of the scattering matrix 
f^KK'i one should calculate either its distribution function or its moments, using the dis- 
tribution d^) of the positional fluctuations. We chose the second procedure, so that the 
magnitude of interest will be the mean square deviation, defined as 

<Ak' = {^KK') - {^kk) , (5) 

where brackets denote the average over the joint probability distribution Yin P {in)) with 
P{in) of the form @. 

In what follows, we assume that the standard deviation a is small. Accounting then for 
that 

-3/2 



1 1 



\n — m + ^n~im\^ \n — m^^ 



{n — my {n — m)^ 

one can expand 5Jnm in Eq. (|45D in Taylor series up to fourth order with respect to 

[n — my 

we are interested in the contribution to cr^j^/ up to the same order. The corresponding 
expressions for /^kk' and A|^^, read 

J ^i(Kn-K'm) / g -^^ gj^g x 

= iV £ |n-m|3 [~2^-- + ^^nm - YQ^nm + Y^^Lj , (6a) 



J X 2 ^i{Kn-K'm) ^-i{Kq-K'p^ 
pq 



X ("-X X - + i^X X^ + ^X^ X^ ^ (6h) 

Carrying out the average in Eqs. (|^), we will collect all terms up to fourth order in a. The 
calculations are rather tedious but straightforward, so we only quote the final results: 
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(7a) 



9J'a 



2^2 



N 



+ 



N 



p,{K) + p;{K') p;{K) + Pj{K') 



+ 



72JV 

N 



Q^{K) + Q^{K')\ + 2[gio(0) + Qi,{K + K' 



where the functions Qy{K) and Pu{K) are given by 



me 



iKm 



\m\ 



Qu{K) = y: 



and the term m = is excluded in the summations. 
For the magnitude of interest, one obtains 



9JV2 



a 



KK' 



N 



p,{K) + p;{K') 



2 

+ 



N 



(7b) 



(8) 



p^{K) + p;{K') p;{K) + p^{K' 



+ 



72JV 

N 



Q^{K) + Q^{K')\ + 2[Qio(0) + Qio(ir + K' 



(9) 



We stress that only terms up to fourth order in a are kept. 

To treat the optical properties of linear aggregates, the important wavenumbers K and 
K' are the smaller ones, namely \K\^ \K\' <^ 1. In this long wavelength limit, the functions 
Pu{K) and Q„{K) can be approximated by their values within the nearest-neighbor frame- 
work, namely Pv{K) ^ 2iK and Qv{K) ^ 2, accounting in the sums (^ only for leading 
terms with m = ±1. Equation (0) then reduces to 



a 



KK' 



36J^ 
N 



{K-K'Ya\l + AQa') + 



1728J2 



N 



-a 



(10) 



Now it can be clearly seen why one should keep the terms up to fourth order in a. The first 
term, being of second order with respect to cr, is equal to zero for the diagonal elements of 
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axK' [note that it is true independently of the magnitude of K since Pu{K) + P^{K) = 0]. 
Thus, for aj^^ one gets 

2 1728 4 , . 

^KK = -^^^ ■ (11) 

This resuh differs from that for uncorrelated diagonal disorder, where the corresponding 
magnitude scales as ci^/A^, being the variance of the site energy distribution [HI]. The 
same behavior appears as well when one simulates off-diagonal disorder by uncorrelated 
randomness in the nearest-neighbor hopping integrals (here, stands for the variance 
of the corresponding distribution). 

With respect to the off-diagonal elements of o"!-^,, we should note that, in spite of the 
fact that the first term scales as cr^, it has an additional suppression factor proportional to 
{K — K'Y ~ A^~^, and thus it may be smaller than the fourth order one. When the first term 
in Eq. ( p!OD dominates crlj^K' cr'^/N^, while in the opposite case one has <Jkx' cr^/N. Both 
results also differ from the scaling law a'^/N found for the other types of disorder [^§,0. 

The origin of the difference discovered lies in that the terms linear in ^„ in the fluctuations 
of the dipolar coupling of site n to the adjacent ones have the same magnitude but opposite 
signs, thus appearing to be correlated, notwithstanding the fact that the fluctuations of 
molecular positions are completely uncorrelated. Due to this feature, they almost (or exactly 
a.t K = K') cancel each other in (^ when summing over n in the long wavelength limit 
{\K\AK'\ 

As it was noted in the previous Section, for perturbative magnitudes of disorder, the main 
effect of (Jkk' is the broadening of the exciton levels. The value of cfrk given by Eq. (pJ]) 
gives the half width of the i^-th exciton state and will serve for determining the latter 
provided that ctxk' remains smaller than the corresponding energy differences \Ek — Ek'\- 
Since the minimum energy difference in the exciton spectrum is between the state with 
K = Q and the next one with K = 2Tt/N , the equality 

(^K=2-k/N, K'=0 = \EK=2Tr/N — Ek'=o\ (12) 
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determines a value of a which separates the ranges of perturbative and nonperturbative 
magnitudes of disorder. In (jl^), Ek is given by Eq. (^). 



IV. DISCUSSION OF THE NUMERICAL RESULTS 

As it was aheady mentioned in the Introduction, numerical simulations of optical prop- 
erties of linear molecular aggregates with off-diagonal disorder generated by Gaussian un- 
correlated fluctuations in the molecular positions yielded different behaviors of the optical 
observables as compared to those for diagonal disorder In this Section, it will be 

shown that the peculiarities found in Refs. can be qualitatively explained from the 

viewpoint of the modified motional narrowing formula (p!oD . In particular, we will focus on 
the dependence of the absorption band width a* and the radiative rate enhancement factor 
on the degree of disorder a. The radiative rate enhancement factor is proportional to the 
number of coherently bound molecules N* , while a* can be estimated from okk replacing 

by A^* |I|,pn]]. Thus, the magic number A^* is, in fact, the unique quantity determining 
the observables we are interested in and, correspondingly, our main goal is to obtain the 



dependence of A^* on a. To do that, we will follow a simple rule earlier proposed in Ref. ||T2 
(see also |ll|]) that works surprisingly well for explaining the corresponding data for diag- 
onal disorder |Tl|,[l2|. This rule simply consists of applying the formula (|12|) to the typical 



localization segment of size A^*. In our estimates of A^*, we will keep only the second term 
in Eq. ([T0| ) since, as will be shown later, it is the major contribution for the parameters used 
in Refs. ||,[10|. 



As a first step, let us take the exciton energy spectrum in the nearest-neighbor approxi- 
mation: Ek = 2 J cos K ^ 2J— JK^. After these simplifications, we arrive at the following 
formula for the number of coherently bound molecules 

^4 \ 1/3 

108 



N*=l—] ^ ^-4/3 (13) 



Estimating now the absorption band width a* as 2a kk by replacing A^ by A^*, we find 
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a* = 2v/l728 Ja^^^ ^ 83 Ja^/^ . (14) 



Recall that for diagonal disorder the respective quantities scale as cr"^/^ 01131 ^^"^ 
0"^/^ P,pi]|-P^, whenever the nearest-neighbor approximation is adopted. It is worth noting 
that the new scaling laws (|T^) and (|T^) follow from the latters (cr^^/^ and o"^/^) simply re- 
placing a hj a^, reflecting the fact that in our case the effective disorder scales as o"^ instead 



of a [see Eq. ([TT|)]. 

The authors of Refs. P,|lO| did not restrict themselves to the nearest-neighbor approxima- 
tion but took into account all dipolar couplings. Using the parameterization ccr", they found 
that their numerical data for the factor of radiative rate enhancement and a* were fitted 
by the sets c = 0.20, a = —1.64 and c = 425J, a = 2.84, respectively. The corresponding 
exponents for the case of diagonal disorder were found to be —0.74 and 1.34, respectively. 
Here, one also observes approximately the two-times increase of the exponents when passing 
from diagonal to off-diagonal disorder. This feature unambiguously shows that the effective 
degree of disorder in the last case scales as a^, what perfectly correlates with our finding 
given by Eq. (|n|). 

Comparison of the numerical fits of Refs. with our results ([TB| ) and (0), obtained 

under the assumption of nearest-neighbor coupling, shows that the numerical cr-scaling of 
both quantities is reproduced reasonably well by the theoretical estimates, better for the 
absorption band width and worse for the number of coherently bound molecules. Similar 
peculiarities are present in the case of diagonal disorder (compare 4/3 with 1.34 and —2/3 
with —0.74). As it was firstly mentioned in Refs. p|JTO[|, including all dipolar couplings 
affects largely the factor of radiative rate enhancement, rising it by more than a factor 2 as 
compared to that calculated in the nearest-neighbor approximation. The explanation was 



done in Ref. [|I6[ and was based on the exact exciton energy spectrum close to the bottom 
of the band 

Ek = -2 Je(3) + JK' (^^ -InK^ , (15) 
where ^(3) = Z^J^i = 1.202. Notice that the presence of a logarithmic term in this 
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equation results in a larger energy level separation as compared to the case of the nearest- 
neighbor model. Thus, a smaller number of the exciton states will be effectively mixed by 
disorder giving finally rise to an increase of A^*. Indeed, if we now use Eq. (^) for finding 
A^* we then arrive at the equation 

lniV*-ln(27r) + 3/2 ~ eTs^ ' 
For the nearest-neighbor model, one should substitute the denominator in the left hand 

side of Eq. (|16D by unity, getting then the previous expression for A^*, Eq. (|T3l). At a fixed 
magnitude of a, keeping the denominator will increase the value of A^*. 

It should be stressed that, rigorously speaking, from Eq. (p!6D does not follow a power- 
like behavior of A^* against a. Therefore, we tried to fit the numerical data of Refs. p|JTO 



relative to the factor of radiative rate enhancement and the absorption band width using 
the following parameterization 

N*\/W b 



In A^* — a 0"^ 



(17a) 



a* = c^. (17b) 

The fits were reached provided a = 0.004, b = 0.12 and c = 148. Note that the numerical 
factor in ( p, 7b|) (c = 148) does not differ drastically from the theoretical value 48v^ ~ 83. 

Finally, we would like to show that for disorder in the interval 0.025 < a < 0.08, used 
in the numerical simulations 0,|l3]; the first term in Eq. ([TOD, neglected by us when making 
estimates, is smaller than the second one. Their ratio at i^' = and K = 2tt/N is equal 
to 7r^(l + 40(7^)/(12A^^(T^). Substituting here A^ by A^* ^ a"^/^, we indeed get for this ratio 
a small mag nitude 7r2(l + ma'^)a'^/^ /12 < 1 when a ranges within the above mentioned 
interval. Thus, our assumption is self-consistent. 



V. CONCLUSION 



The motional narrowing effect in one- dimensional Frenkel chains with off-diagonal dis- 
order arising from Gaussian fluctuations in the molecular positions is found to be different 
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from that for diagonal disorder and for off-diagonal disorder with uncorrelated randomness 
in the nearest-neighbor couplings. Such distinction is due to the fact that the fluctuations 
of the dipolar coupling of a given molecule to the adjacent ones are correlated, even if the 
fluctuations of the molecular positions are completely uncorrelated. Thus, this type of dis- 
order cannot be modeled by uncorrelated randomness in the nearest-neighbor interactions. 
The estimates of scaling of the optical observables with the degree of disorder based on the 
new motional narrowing law are found to be in qualitative agreement with those obtained 
previously in numerical simulations. 
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